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Abstract 



We investigate the electronic structure of the helium atom in a magnetic field between B = and 
B = 100a. «.. The atom is treated as a nonrelativistic system with two interacting electrons and a fixed 
nucleus. Scaling laws are provided connecting the fixed-nucleus Hamiltonian to the one for the case 
of finite nuclear mass. Respecting the symmetries of the electronic Hamiltonian in the presence of a 
magnetic field, we represent this Hamiltonian as a matrix with respect to a two-particle basis composed 
of one-particle states of a Gaussian basis set. The corresponding generalized eigenvalue problem is solved 

^^ . numerically, providing in the present paper results for vanishing magnetic quantum number M = and 

even or odd z-parity, each for both singlet and triplet spin symmetry. Total electronic energies of the 
ground state and the first few excitations in each subspace as well as their one-electron ionization energies 
are presented as a function of the magnetic field, and their behaviour is discussed. Energy values for 
electromagnetic transitions within the M — subspace are shown, and a complete table of wavelengths 

Qs ' at all the detected stationary points with respect to their field dependence is given, thereby providing a 



basis for a comparison with observed absorption spectra of magnetic white dwarfs. 



1 Introduction 



Since the astrophysical discovery of strong magnetic fields on the surfaces of white dwarfs (< 10 5 Tesla) and 
_Cd \ neutron stars (~ 10 8 Tesla) the interest in the behaviour and the properties of matter in strong magnetic 
fields has increased enormously. In the case of atoms in strong magnetic fields most of the literature is 
concerned with the hydrogen atom. The eigenvalues and the eigenfunctions of the hydrogen atom are, 
therefore, known very precisely and for many excited states |l]H7]]. For the hydrogen atom in a strong 
magnetic field it was possible to perform a comparison of the theoretical predictions resulting from ab initio 
computations with the data obtained from astronomical observation. The results provided an overwhelming 
evidence for the existence of hydrogen in the atmosphere of the corresponding astrophysical objects ||. 
However, there are several spectral features and structures which can not be explained by hydrogenic spectra 
like for example in the spectrum of the magnetic white dwarf GD229 [p}-|l0||, which leads to the conjecture 
that there are further components in the atmospheres, i.e. atoms with more than one electron. 

* Permanent address: Department of Physics, Univ. of Athens, GR- 15771 Athens, Greece 



So far our knowledge about atoms with more than one electron in strong magnetic fields is very sparse. 
Most of the ab initio computations on more electron systems deal with two-electron systems, i.e. the helium 
atom, the hydrogen anion and the helium- like kations K, 11 14J . For two-electron atoms the investigations 



cover much smaller parts of the spectrum than it is the case for the hydrogen atom, and their accuracy is 
considerably poorer. To perform a comparison with astrophysical observation, however, accurate transition 
energies for a large set of field strengths have to be available which requires even more accurate data for the 
total energies. In particular in the relevant regime of intermediate field strengths, for which the diamagnetic 
energies and the Coulomb energies are of the same order of magnitude, there does not exist sufficiently 
accurate data for a large number of excited states which would allow this comparison. The reason is that for 
example the numerical basis set methods established so far have their starting point either in the low-field 
regime or in the high-field regime and are, therefore, especially adapted to their corresponding regimes but 
fail to be effective in the opposite regime, thereby leaving a gap in a certain regime where none of them 
provides a good convergence. 

The scope of the present paper is to introduce one uniform basis set method which is capable of accurately 
describing two electron systems for arbitrary field strengths. We present numerical results for the energy 
levels of the ground state and a considerable number of excited states of the helium atom in a magnetic 
field ranging from B = Oa.u. to B = 100a. u. (B = la.u. corresponds to 2.35 • 10 5 Tesla). 

In order to perform a detailed comparison of our numerical data with already existing ones, let us first 
provide an overview of the results present in the literature so far. Almost all of the results we mention here 
are variational, i.e. they provide upper bounds for the exact energy values. However, we have to distinguish 
between Hartree-Fock calculations and calculations which take into account the electronic correlation. 

The first Hartree-Fock calculations on helium in a strong magnetic field have been performed in 1976 
by Virtamo [|1|. Virtamo has provided the global ground-state energies of helium in very strong fields 
(B = 80a. u. up to B = 8 x 10 4 a.u.). In this regime the lowest energy is achieved by aligning both spins 
antiparallel to the magnetic field which means that the global ground state in question is a spin triplet 
state. The same state is also investigated in the work of Proschel et al.[0]. They use a Slater-determinantal 
approach starting from Landau-levels in order to cover a slightly broader regime of field strengths as Virtamo 
(B ~ 2. la.u. to B ~ 2.1 x 10 4 a.u.). Hartree-Fock results for several other states besides the ground state 



of helium in the high-field regime are provided in the works of Ivanov [14, 16 where the singlet and triplet 
states are considered for positive total z-parity and magnetic quantum numbers m = 0, —1 as well as for 
negative z-parity and m = 0. Here the magnetic field ranges from B = Oa.u. to B = 100a. u. occupying 11 
field strengths and thus includes both the low-field and high-field regime as well as the intermediate regime. 
On an even finer grid of field strengths (34 values between B = 8 x 10~ 4 a.u. and B = 8 x 10 3 a.u.) the 
energies of the triplet ground state and several triplet excited states but no singlet states are given in the 
work of Thurner et al. p3| . In the latter work there exists a gap in the list of eigenenergies for several states 
in the vicinity of B ~ la.u. because the ansatz of the applied approach changes in the intermediate regime 
from a spherical to a cylindrical symmetry and thus fails to provide accurate results in the intermediate 
regime. An important work is also the Hartree-Fock study of Jones, Ortiz and Ceperley p7| . They present 
the HF-energies for the global low-field ground state for various field strengths between fi = 8x 10 _4 a.u. 
and B = 8a. u. and thus also address the intermediate regime. They additionally provide energies for several 
excited states (but those results are crude approximations providing no upper bounds). 

In comparison to the considerable number of Hartree-Fock investigations there are much less data avail- 
able on fully correlated calculations of helium so far. Mueller, Rau and Spruch [18| obtain variational 



upper-bound estimates of the binding energies of the triplet states with negative z-parity and the magnetic 
quantum numbers m = and m = — 1 and of the singlet ground state which has positive z-parity and 
m = 0. Their field strengths range from B ~ 4.2a.u. to B ~ 2.1 x 10 4 a.u.. A similar high-field regime 



{B = 4. 2a. u., 42a. u., 420a. u.) is addressed by the work of Vincke and Baye |T^] which also presents varia- 
tional estimates for binding energies. They consider the singlet and triplet states for positive z-parity and 
m = 0, — 1, —2. The first correlated calculations in the intermediate regime have been provided by Larsen 
[ pc[| who has given the energies of the singlet states with positive z-parity and m = 0, —1 or with negative 
z-parity and m = as well as the triplet ground state for the four field strengths B = 0.2,0.5, 1.0, 2.0a.u.. 
Park and Starace pi]] have computed upper and lower bounds for the singlet ground state in the low-field 
regime between B = Oa.u. and B = 0.15a.u.. Recently, Jones, Ortiz and Ceperley |22| have applied a 
released-phase quantum Monte Carlo method to the helium atom. Energies for spin triplet states with both 
positive and negative z-parity and m = 0, — 1 are given for field strengths from B = Oa.u. up to B = 8. Oa.u.. 
Though in principle the results of those investigations are variational their numerical values must be handled 
with care when using them as upper bounds for the helium energies. The reason is that due to the statistical 
character of their method the results possess error bars which might trespass the exact values of the helium 
energies. Very recently the finite-element technique has as well been used to calculate energies for several 



singlet and triplet states of helium in a magnetic field [23]. Also this approach does not provide a significant 
improvement. 

In the present work we use a Gaussian basis set method which permits to perform fully correlated 
calculations on singlet and triplet states of helium with arbitrary spatial symmetry. This basis set method 



has been developped by Schmelcher and Cederbaum [24] for molecules and has already successfully been 



applied to the hydrogen molecule ion [25] and to the neutral hydrogen molecule |26|,|27j in strong magnetic 
fields. We translate this basis set in order to be applicable to ab initio calculations of atoms in strong 
magnetic fields. We remark that using for example a Hylleraas basis set yields more accurate results in the 
field-free or weak-field case due to explicitely correlated orbitals. However, it is not efficient and accurate in 
the presence of a strong magnetic field due to the spherical symmetry of the exponentials [28]. In contrast 



to this our basis set described in detail in sec. Ill below can be adapted to any field strength. 

All the helium states considered are classified according to a maximal set of conserved quantities, which 
we choose to be the total spin S 2 , the z-component S z of the total spin, the total spatial magnetic quantum 
number M and the total spatial z-parity Ii z . For twenty field strengths from B = Oa.u. up to B = 100a. u. 
numerical helium energies for singlet and triplet states are given for M = and for U z = 0, 1. In each of 
these two subspaces we present the energies of the ground state and the first five excited states for singlet 
and four excited states for triplet spin symmetry. The accuracy of our method in the field-free case can 
be determined by comparison of our results with the very precise field-free values calculated by Braun et 
al.|2i|, Accad et al.[3C] or by Drake et al.[31], and our relative deviation ranges from 10 -6 to 10~ 4 . We do 



not observe this accuracy to drop considerably with increasing magnetic field, which means that our basis 
set method produces accurate results in particular in the intermediate regime. Both with respect to the 
achieved relative accuracy as well as with respect to the number of excited states our investigation provides 
therefore valuable results. Data for M/0 will be presented in a future work. The stationary components 
of the transitions obtained in the present investigation provide an important part of the data which have 
very recently been used to show the presence of helium in the atmosphere of the magnetic white dwarf 
GD2290. 

In detail we proceed as follows. In section II we discuss the influence of the finite nuclear mass, i.e. 
the scaling relations which connect the Hamiltonian for a finite nuclear mass to a Hamiltonian with infinite 
nuclear mass. In all of the remaining sections, we consider the nonrelativistic case with infinite nuclear mass 
and an electron spin g- factor equal 2. Section III starts with the nonrelativistic Hamiltonian of a helium 
atom with infinite nuclear mass in a magnetic field. The corresponding symmetries of this Hamiltonian are 
discussed and serve for constructing a suitable basis set. The matrix representation of the Hamiltonian with 
respect to this basis gives rise to an eigenvalue problem. The results of its diagonalization are discussed 



in section IV where we present total energies, ionization energies and transition energies as well as the 
wavelengths of their stationary components, all of them given under the assumption of infinite nuclear mass. 

2 Finite nuclear mass scaling relations and further corrections 

We use a nonrelativistic approach to the helium atom. This is justified because the relative changes in 
the helium energies due to relativistic effects are smaller than the relative accuracy of our nonrelativistic 
energies |32]. Additionally for simplicity we use an electron spin g-factor equal 2 throughout the paper. The 
following scaling laws as well as all of our numerical results can trivially be adapted to a g-factor of any 
desired accuracy by multiplying every occuring spin operator or eigenvalue with |. 

Though our electronic structure calculations will be carried out with the assumption of an infinite nuclear 
mass, they can, by means of a scaling law, be translated into results with a finite nuclear mass. The latter 
are necessary in order to achieve the accuracy for a detailed comparison with astrophysical data. The idea 
to introduce such a scaling relation is not new but has already been applied to one-electron systems earlier 
[|3q] . In order to be demonstrate how this idea is applied on the helium atom, we will start with the full 
Hamiltonian describing a neutral system with two interacting electrons with masses 1 and charge —1 in 
the Coulomb field of a nucleus with charge two and finite mass Mq in a magnetic field B. Assuming a 
vanishing pseudomomentum of the center of mass (see refs. [|34|-36 for a pseudoseparation of the center 



of mass motion for neutral systems in a magnetic field) one obtains the following exact pseudoseparated 
Hamiltonian (which is established in the symmetric gauge A(r) = 2 B x r) 

2 ' 1 o 1 . 1 , 2 \ 1 



He = E U^ + ^ B - 1 * + ^( BXr <) 2 -O + B - S + 



2(jl ' 2/i' 8/j, ' |rj| / |r 2 -ri[ 

O- E (piPj ~ P*( B x t j) + j( B x r *X B x r i)) > C 1 ) 



2M . , 

where \i = M ° 1 is the reduced mass and // = M °^ 1 (we use atomic units throughout the paper). We 
take into account the dominant part of the finite mass corrections by introducing these quantities //, // and 
neglect the mass polarization terms represented by the sum over (i ^ j). Our basic idea is now find a scaling 
law joining the spectrum of the remaining part H(Mq, B) in the first line of the Hamiltonian H e with the 
spectrum of the infinite-mass Hamiltonian H(oo, B) at a suitable different field strength B. The neglect of 
the mass polarization terms of the sum over [i ^ j) in eq.(||) is justified by the fact that they are expected 
to provide a smaller correction to the total energy than the corresponding diagonal mass corrections, i.e. 
the replacement of the masses by reduced ones. 

We will show how it is possible to establish a relation between the operators H(Mq,B) and H(oo,B) 
themselves which is even more fundamental than only a relation between the corresponding spectra. The 
relation in question cannot be a simple Unitarian because such a transformation would leave the spectrum 
invariant. Instead, we will show that there is an unique way how to represent H(Mq,B) in the form 
UH(Mq, B)U~ 1 = aH(oo,B) + 0, where a is a number and is a suitable operator commuting with the 
Hamiltonian. 

Before factoring out a suitable prefactor a the terms of H(Mq,B) must be transformed in such a way 
that their mutual ratios possess the correct infinite-mass values, i.e. all ratios must neither contain /x nor //. 
The misproportion between the kinetic energy and the Coulomb energy terms is the only one which cannot 
be absorbed in the parameter B or be removed by separating a suitable additive operator (3 commuting with 
the Hamiltonian. Thus we are obliged to introduce the canonical scale transformation r — > -r, p — > up of 



the coordinates themselves. We realize this transformation by the Unitarian U = e ' - 



UH{M ,B)U~ l = fi 



E f k + ^B x r, 
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^(XP + PX^ y ielding 
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2/i'' 



-B • L r + B s, 



(2) 



The operator Sj of the spin degree of freedom is not affected by the above transformation, and L. = rj x pj is 
also unchanged because the scaling factors of r-j and p« cancel each other. The expression in square brackets 
is already an essential part of the desired infinite-mass Hamiltonian at the adjusted field strength B = —?B: 
only the Zeeman orbital and spin term at the field strength B are missing. These operators can be provided 
by hand, which must be repaired by subtracting the same operators at the end where they almost cancel 
exactly the original Zeeman and spin terms leaving only a contribution of order jj^ : 

2 



UH{M ,B)U 
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E(^ + ^B.l l + i(Bxr l ) 2 -^ + B 
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|r 2 -ri 



E 



-z^B • 1* 



fjB ■ Si + — B • U + B 
2/i' 



(3) 



Therefore we obtain 



UH{M 0l B)U- 1 = /i-7Y(oo,i?//x 2 )-— B-^(1 J + Sj ) (4) 



In particular, the spectrum of H(Mq,B) is identical with the spectrum of the unitarily equivalent operator 
represented by the right side of (0). The latter spectrum can be simply connected to the spectrum of 
H(oo, B//j?) itself because the operators ^ B • lj and ^ B • Sj commute with the Hamiltonian: we will use 
in the following section an angular momentum and spin basis, in the case of which the additional operator 
— jg-B • J20-i + s i) gives rise to a trivial energy shift of order jj^. 



3 Symmetries, Hamiltonian and Basis sets 

3.1 Symmetries and Hamiltonian 

In the following we assume the magnetic field to point in the +z-direction. Then the Hamiltonian H(oo, B) 
for infinite nuclear mass at given field strength B (in the following we will omit the argument 'oo') reads 
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H = E(5P? + ^« + 



i=i 



[ ixi + yf) 



W\+Bs z 



+ 



|r 2 -ri| 



(5) 



The sum contains the one-particle operators, i.e. the Coulomb potential energies — "tjtt of the electrons in 
the field of the nucleus as well as their kinetic energies, already splitted into the part 2P 2 > the Zeeman 



term ^Blzi and the diamagnetic term 



{xl + Vi) 



and their spin energies Bs zi . The two-particle operator 



ra—ri 



represents the electron-electron repulsion energy. 



There exist four independent commuting conserved quantities: the total spin S 2 , the z-component S z 
of the total spin, the z- component L z of the total angular momentum and the total spatial z-parity Ii z . 
In the following calculations we consider separately each subspace of a specified symmetry, i.e. with given 
eigenvalues of S 2 , S z , L z and Il z . 



3.2 The underlying one-particle basis set 

The key ingredient of our basis set method is the Gaussian one-particle basis set which is our starting point 
for the construction of spatial two-particle states. According to the azimuthal symmetry with respect to the 
magnetic field axis, cylindrical coordinates are suitable for representing the one-particle basis functions 

$i(p,<p,z) = p n p,z n ^e-^ p2 -^ z2 e im ^ i = l,...,n , (6) 

where on and Pi are positive nonlinear variational parameters and the exponents n pi and n zi obey the 
following restrictions: 

n Pi = \mi\ + 2ki ; ki = 0, 1,2, ... with mi = ... — 2, — 1,0, 1,2, ... (7) 

n zi = n zi +2k; h = 0, 1,2, ... with ir zi = 0, 1 (8) 

The basis function <3?j is an eigenfunction of the z-component of the angular momentum with eigenvalue rrn 
and an eigenfunction of the z-parity with eigenvalue (— l) 71 " 2 *. The Gaussian-like expression p\ m i\e~ aip is 
identical with the p-dependence of the lowest Landau-state in the field B if we choose Qj to be -j and thus 



represents an adjustment to the existence of the magnetic field (see e.g. eq. (8) in |36|]). The monomials 
p2ki an( j z n zi are su ^ ame f or describing excitations. The flexibility of our basis set which permits us to 
choose suitably the values of the nonlinear parameters a« and Pi is one of its major advantages: For low field 
strengths an isotropic choice a» = Pi will be reasonable where the a^, Pi cover a regime which allows the $j 
to optimally approximate Slater- type orbitals. At high field strengths, however, where the magnetic field 
destroys the spherical symmetry of the problem, an isotropic basis set method would be inefficient. Here we 
choose the distribution of the a^ to be peaked around -j whereas the Pi are well-tempered in a large regime. 
The best choice of the «j and Pi has been computed by the requirement to solve optimally the one- 
particle problem of the if -atom or the He + -ion in a magnetic field of given strength B. To achieve this, 
we applied the following optimization procedure. For any given one-particle subspace (m,ir z ) which will 
be involved in the later on two-particle configurations we have chosen a suitable number, typically 20, of 
functions of the type @ with the same one-particle quantum numbers (rrij = m, ir zi = ir z ) but different 
starting values Oj and Pi. With respect to this basis we built up the overlap matrix and the matrix of the 
one-particle Hamiltonian, representing a generalized eigenvalue problem for the one-particle energies. Now, 
we determined systematically the values of the «j and Pi which minimized the energy eigenvalues of the one- 
particle ground state or a desired excited state. Our tool was a repeated reconstruction and diagonalization 
of the matrices directed by the pattern-search- method applied in the {aj, Pi} space. We remark that due to 
the large number of parameters it is very time consuming to find the best selection of parameters {/c,,/,}. 
Additionally the starting values for the nonlinear parameters «j and Pi have to be chosen carefully in order 
to find a "good" minimum close to the global one on the complicated hypersurface. 

3.3 The two-particle basis set 

Our basic idea for solving the time-independent Schrodinger equation for the Hamiltonian (||) is to construct a 
basis set of suitable two-particle states \q) each of which is itself an eigenstate of the conserved quantities with 
eigenvalues M,U Z , S, S z , respectively. With respect to this not necessarily orthonormal basis we construct 
a matrix representation H pq of the Hamiltonian and define an overlap matrix S pq by 

H pq := {p\H\q) (9) 

S pq := (p\q) (10) 



where the states \p), \q) lie in the same subspace of given eigenvalues of M, IL Z , S, S z . By construction, the 
matrix H_ is hermitian and the matrix S is hermitian and positive definite. Additionally to these general 
properties, in the special case of the basis sets we use, all the matrix elements turn out to be real. By solving 
the finite-dimensional generalized real-symmetric eigenvalue problem 

(R-ESJ-c = (11) 

we obtain eigenvectors c whose corresponding eigenvalues E are strict upper bounds to the exact eigenvalues 
of the Hamiltionian (g) within the given subspace. 

We choose our two-particle basis functions to be direct products of a pure spatial part \ip q ) with a pure 
spin part \Xq)' 

\q) = h^)®lx«) (12) 

The spin part is assumed to be one of the usual orthonormal singlet or triplet spin eigenstates. This means 
that the overlap matrix is nontrivial purely due the spatial part, 

S pq = (i> P \%) , (13) 

and each matrix element of the Hamiltonian decomposes in a spatial and a spin contribution: 

H pq = {tpp\ H spat \lpq) + Spq {Xp\ H spin \\q) (14) 

The spin contribution (Xp\Hspin\Xq) is trivial and equals +B ■ S z , according to (|5|), and therefore it vanishes 
for singlet states or is either —B,0 or -\-B for triplet states. Due to its prefactor S pq it can simply be 
absorbed as a shift in the energy eigenvalue E in eq.([ll|). Since a spin singlet state is antisymmetric and 
a triplet state is symmetric with respect to particle exchange, the spatial part associated to a spin singlet 
state must be symmetric whereas the spatial part associated to a spin triplet state must be antisymmetric. 
In contrast to the spin part the spatial part is by no means trivial and it is an art to find a suitable 
finite-dimensional basis set for accurately approximating the Hilbert space of bound states of H spa t. This 
is the reason why we spent so much effort on providing a powerful one-particle basis set (see sec. &%) ■ It 
serves now as a good starting point to compose the desired two-particle states. We choose a two-particle 
basis state to be 

|V 9 > := &I&t|0> i = l,...,n , j = i,...,n , (15) 

where b\ is a creation operator of the i-th one-particle basis state \i) = b\ |0) whose position representation 
is given by (|6|). This means that we will treat the helium atom by a full Configuration Interaction (full 
CI) approach within the one-particle basis set of states \i). Depending on whether the spin part \xq) is a 
singlet or triplet state, the operators b\ must be bosonic or fermionic, respectively. 

Now, in order to establish a basis set of two-particle states spanning the subspace with the given pure val- 
ues of the total magnetic quantum number M and the total z-parity II^, we must select all the combinations 
i,j with 

rm + m j = M (16) 

mod(TT zi + TT zj ,2) = Ii z , (17) 

yielding a dimension N of the constructed two-particle basis set which is in general smaller than — ^ — -. 
We remark that for triplet spin symmetry those states in ( |l5| ) with i = j fail to exist which means that the 
dimension N for triplet subspaces is in general smaller than for singlet subspaces. 



3.4 Matrix elements 

In order to calculate the matrix elements of the Hamiltonian (||) , we must rewrite its spatial part in second 
quantization, H spa t = Hi + Hjj, where Hj and Hjj denote the second-quantized counterparts of the familiar 
one- and two-particle operators whose position representations read 

Hi(p,r) = ip 2 + iB-l + |(x 2 + y 2 )-^ (18) 

fljr(ri,r 2 ) = ]^hr\ (19) 

The next step now is to calculate the spatial matrix elements according to (|l4|). With 1^} := b\bj |0) and 
\ip p ) := b\b\ |0) a straightforward calculation leads to 

(^|^> = (i\k) 010 ± (i\l) (j\k) (20) 

> p |flj|^) = (i\Hi\k) 010 ± (%\Hi\l) (j\k) 

+ (j\Hj\l) (i\k) ± (j\Hj\k) (i\l) (21) 

> p |^/r|^) = (%j\H n \kl) ± (ij\H„\lk) , (22) 



where |ij) := |i) <8> |j) and where the sign '±' stands for '+' in the singlet case and for '— ' in the triplet 
case. 

All the — g — different one-particle matrix elements (i\Hj\k) are relatively easily evaluated (see Appendix 
^). Each matrix elements (i\Hj\k) equals a prefactor depending on the parameters of \i) and \k) times the 
one-particle overlap (i\k), reflecting the fact that the one-particle operators do not only conserve the total 
magnetic quantum number M and the total z-parity II^ but also the magnetic quantum numbers and 
z-parities of the individual one-particle states. We remark that the mentioned prefactors in the matrix 
elements of the Zeeman term do not even depend on any basis state, rendering them purely proportional to 
the overlap, (i|;jB • l\k) = 2 rn kB (i\k). Thus the Zeeman term just gives rise to a shift in the energies. 

The only coupling between two-particle states composed of one-particle states with different combinations 
of magnetic quantum numbers or z-parities arises due to the two-particle operator of the electron-electron 
interaction. In contrast to the one-particle matrix elements the evaluation of the matrix elements (ij\Hjj\kl) 
is by no means trivial (see Appendix |C]). An additional problem is the large number of different two-particle 
matrix elements which is of the order — *-= — - rather than —^ — • This means that a sophisticated and 
detailed analysis of the analytical representation of the two-particle matrix elements by means of series of 
different representations of hyper geometric functions has been necessary in order to achieve an effective 
numerical implemetation and an acceptable CPU time for performing a calculation for one given field 
strength. Details on the representation of the quantities (ij\Hjj\kl) are given in Appendix |Q. 

4 Results and discussion 

4.1 Spectroscopic notation and properties 

Before presenting the numerical results of our calculations we shall explain our spectroscopic notation in the 
presence of the field as well as its correspondence to the field-free notation. According to the four conserved 
quantities M, H z , S 2 and S z we denote a state by v s + M^~ l > z where (25 + 1) is the spin multiplicity and 
v = 1, 2, 3, ... is the degree of excitation within a given subspace. If obvious, we will omit the index S z in the 
following. In the present paper we will investigate and show results for the subspaces 1 + , 3 + , 1 _ , 3 CP. 



For vanishing field, there exists a one-to-one correspondence between our field notation and the common 
field- free notation n s + Lm (see Table 1). It is given by the fact that for bound helium eigenstates of L 2 and 
L z we have (— l) nz = (— 1) L+M . The latter relation is not in general true for multi-electron systems because 
for a configuration containing two angular momenta l\,l2 we have (— l) nz = (— l)'i+fe+-kf_ For helium and 
B = 0, however, all doubly excited states lie in the continuum |37j which means that all the helium states 
below the one-particle ionization threshold T(B = 0) = — 2.0a.u. must be one-particle excitations with 
l\ = and l-i = L. 

We emphasize that the mentioned one-to-one correspondence for B = does not contradict the fact that 
the L 2 -symmetry is higher than the Il^-symmetry: At a finite field strength the former is broken whereas the 
latter is still valid. The field free notation becomes meaningless for finite field and could only be maintained 
as a labelling device for the energy eigenstates because a given v s + M^" 1 ' z state at a finite field develops 

in a unique way from the equally labelled v s M^~ l > z state at B = which, in turn, is identical with one 

unique n s Lm state. We remark that energy curves within a subspace of given symmetries s M' -1 ' 



are not expected to cross [38|, whereas crossings between curves of different subspaces are allowed. 



In Table 1 we also provide the energies of the field free states which have been very precisely calculated 



by Accad et al. |3(| or by Drake et al. |?T| . We want to keep track of the energetical order of the states for 
several reasons. The first reason is that we need the energies for associating the energy quantum numbers 
v to the states. It would not be sufficient to use the same counting n as in the field free case because there 
exist different states with the same n but different L which possess the same n z -symmetry, as is evident 
already for the states 3 1 + and 4 1 + . 

The second reason is to point out the approximate degeneracy of the field free states with the same 
energy quantum number n but different L. Whereas the electron-electron interaction is only able to slightly 
perturb the otherwise exact degeneracy of those states, we will see that the magnetic field will completely 
remove this degeneracy. This effect can even be observed for the behaviour of the corresponding states in the 
H atom and thus is primarily an effect of the magnetic field alone rather than of the two-particle character 
of the He atom. It occurs in addition to the well-known removal of the degeneracy of states with the same 
quantum number L but different M whose energies split for finite field in (2L + 1) different energy values. 

The properties of the field free states of helium summarized in Table 1 serve as a good starting point for 
presenting our data for finite fields. 

4.2 Aspects for the selection of basis functions 

In order to obtain accurate results by our basis set method two major difficulties have been overcome by an 
optimal selection of basis functions. The first difficulty is the limited number n of one-particle functions from 
eq.(P) which can be used to describe the exact wave function. The second difficulty, which is not completely 
independent from the previous one, is to describe electronic correlation by using a basis composed of one- 
particle states. 

One manifestation of correlation is the fact that different electrons avoid occupying the same region in 
space which we expect to be the consequence of the electron-electron repulsion as well in a magnetic field. 
In particular, we decribe this by two electrons tending to occupy two regions in opposite directions from the 
nucleus which corresponds to the situation (f2~ fi = tt. The Coulomb interaction , _ r , breaks the indepen- 
dent conservation of the ^-components l z i,lz2 of the two angular momenta, only leaving the sum L z = l z \+l Z 2 
as conserved quantity. A Fourier representation of the angular part of any fully correlated two-particle wave 
function ^(ri, Y2) is given by ^(^i, v?2) = J2m 1 ,m 2 A mim2 e l ( mi ' fil+m2 ' fi2 \ In order to obtain an eigenfunction 
of L z we must demand the constraint m\ + rri2 = M, yielding ^(^1,922) = e tMipi Y, m ^-Af-m^e 1 " 1 ^ 2 "^ 1 ^. 
This expression contains angular correlation as can already be seen explicitely by considering the lowest 



cosine term B m cosm((p2 — (pi) contained in the series above: assuming B m to be negative, ^((fi, ^2) is 
largest for y?2 — <-P\ = 7T and lowest for y?2 — ^1 = 0. Therefore, it appears fruitful to choose one-particle wave 
functions with opposite magnetic quantum numbers ±1 and ±2 or even higher angular momenta combining 
to a total magnetic quantum number M = in order to describe well the angular correlation. Correlation 
can also been described by wave functions possessing a node at the nucleus position and thus allowing the 
two electrons to be located in opposite directions with respect to the nucleus. 

The maximum dimension iV of the Hamiltonian matrix is limited by CPU and storage resources. Through 
a very efficient implementation (see appendices |A|.|B|,|C|) of the matrix elements and by optimized storage 
usage we were able to push the number of two-particle basis functions to N ~ 4300 which represents the 
limit with respect to linear dependencies of similar basis functions resulting in instabilities of the numerical 



diagonalisation of the generalized eigenvalue problem (11). 

Let us discuss the strategy for the selection of basis functions for the example of the subspace + . We 
first focus on the basis functions optimized for a nuclear charge Z = 1, i.e. the hydrogen atom. We selected 
31 functions with the symmetry + , among which 13 have a p-exponent equal 2 instead of in order to 
describe correlation with the aid of their nodes. For further improving the description of correlation we used 
also each 13 functions with m' Kz = ±1 + and m' Kz = ±2 + as well as each 13 functions with m nz = 0~ and 
m? z = ±1~. The first excited state is already rather well described by this basis set but is still improved 
by adding basis functions optimized for excitations of the H atom in a magnetic field. Whereas for the first 
excited state correlation effects are still important the higher excitations are, like in the field free case, more 
and more dominated by one-electron excitations for which automatically the two electrons are spatially 
separated. Therefore, higher excitations were throughout described by adding functions with values ki = 1 
or li = 1 in eqs.(|3) and <JS|) which are subsequently optimized to describe the corresponding higher excitations 
within the + subspace of hydrogen. 

In order to describe higher excitations of the He atom the optimization for Z = 1 is sufficient because 
the nucleus with Z = 2 is screened by the inner electron. For the ground state and the first two excited 
states, however, we have observed important contributions also from functions optimized for Z = 2. Here 
we used a similar selection scheme as for the functions optimized for Z = 1, but only half the number of 
functions. 

Altogether we arrive at a number of n = 244 different one-particle states, from which Aq = 4378 two- 
particle states for singlet and A3 = 4288 two-particle states for triplet spin symmetry can be composed 



according to eq.(|l5|) 



In a similar manner, we built up the basis for the 0~ subspace by involving functions optimized for Z = 1 
as well as for Z = 2. We observe the effect of the latter optimized functions to be less significant in the 0~ 
case than in the + case. This appears natural since in any two-particle state with 0~ symmetry at least one 
electron is in an excited one-particle state and is only attracted by a screened nucleus with effective charge 
closer to Z = 1 than to Z = 2. Due to this reason we used altogether n = 195 one-particle states which is 
slightly less than in the + case. According to eq.(15) we obtain Ai = A3 = 3600 two-particle states for the 



singlet and triplet subspace, respectively. The numbers N± and A3 do not differ due to the fact that the 
identical spatial one-particle quantum numbers which are always forbidden for the triplet two-particle states 
do not occur for the 0~ singlet subspace either because two different one-particle z-parities are required to 
form an odd total z-parity. 

In the following subsection we discuss the results of our He calculations. Their accuracy is estimated by 
comparison with the field free data. If available, we provide also a comparison of our data for finite field 
strengths with the literature. 
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4.3 Energies for finite field strengths 
4.3.1 Results for M = and even z-parity 

a) Results for the singlet states u 1 + 

For the singlet subspace u 1 + we present the ground state and the first four excitations, i.e. 1 < u < 5. 
In the low-field and part of the intermediate regime the 1 1 + state is the global ground state. The numerical 
results for the total energy of the 1 1 + state as a function of the magnetic field are shown in Table 2. We 
observe that apart from B = and B = 0.08 our results for the energies are throughout lower and thereby 
better than the best ones given in the literature. We remark that the 1 1 + state is the state which causes the 
most difficulties for an accurate description by our method. The first reason is that it is the only state for 
which both electrons considerably occupy the one-particle ground state which forces the two electrons to be 
close to each other in a narrow domain of space and which thus gives rise to a relatively strong contribution 
of electronic correlation. The second reason is that the mentioned one-particle ground state, which would 
be a Slater-type orbital for B = 0, possesses a cusp at the origin which is difficult to obtain accurately by 
a superposition of Gaussians. This second effect, however, can be expected to become less important with 
increasing field strength since the cusp in the direction perpendicular to the field axis is smoothened out by 
the increasing dominance of the magnetic field. 

The overall increase of the total energies (see Table 2) with increasing field strength has its origin in the 
strongly increasing kinetic energy in the presence of the external field. The total binding energy of the two 
electrons can be obtained from the total energies E(B) in Table 2 by subtracting the minimal energy B of 
two free electrons in the field, yielding E(B) — B (here and in the following all ionizations are considered 
with fixed quantum numbers). 

However, for an analysis of the electronic structure the one-electron ionization energies for the process 
He — > He + + e~ are much more sensitive. The corresponding one-particle ionization threshold T{B) is 
provided in the fourth column of Table 2. This threshold T can easily be obtained as the sum of the 
lowest Landau-energy -^ of the ionized electron and the total one-particle energy E\ l{Z = 2) of the other 
electron in the Coulomb field of the nucleus with charge Z = 2 in the presence of the magnetic field. 
This quantity E\ l(Z = 2), in turn, can be received from the highly accurate values for the one-particle 
binding energy EL'JZ = 1) computed by Kravchenko et al. ||. First we extract the total energies for 
Z = 1 from {—E hi ^ d {Z = 1)) := E\ l{Z = 1) — -=-, and then we use the nuclear charge scaling relation 
E(Z,B) = Z 2 E{Z = l,B/Z 2 ) in ref.g, yielding 

T(B) = B-4E$ ld (Z = l,B/4) (23) 

This threshold T lies essentially closer to the values E(1 1 + ) than the higher threshold B for two-particle 
ionization. We consider it in fact to be the optimal reference point for the total energies, and as expected 
the total energies of the excitations u 1 + which are given in Table 3 approach closer and closer to the value 
T with increasing excitation. In agreement with the mentioned considerations about the special difficulties 
in achieving accurate results for the ground state we obtain much preciser values for the excitations within 
the same basis set. 

The one-particle ionization energies for all the states v 1 + , 1 < v < 5, i.e. the values of E(B) — T(B), 
are given in Fig.l as a function of the magnetic field. A logarithmic scale on the energy axis is necessary 
for covering the three orders of magnitude for the different states. We observe that though the ground 
state 1 1 + becomes monotonically stronger bound with increasing field strength this is not in general the 
case for the excitations. Below S~0.005a.-u. none of the energies differs considerably from its field free 
value. Between i?~0.005a.u. and B^O.la.u. a rearrangement takes place which is caused by the increasing 
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dominance of the magnetic forces over the Coulomb forces. The first effect to be observed is that the 4 1 + 
state leaves the energetical vicinity of the 3 1 + state. According to Table 1, for B = these states would 
coincide with the states 3 1 Sq and 3 x Dq, respectively, and thus both correspond to the same energy quantum 
number n = 3. The degeneracy of these two states is only slighty disturbed in the field-free case whereas a 
finite field removes this approximate degeneracy completely above 5~0.01a.it.. 

At -B~0.08a.it. the states 4 1 + and 5 1 + experience an avoided crossing. It can be recognized well even 
though the strong curvature of the graphs would make a finer grid of calculated energy values on the field 
axis desirable. The present grid is, nevertheless fine enough to observe that the calculated points are well 
aligned, confirming the good convergence of our calculations. 

b) Results for the triplet states v 3 + 

For the triplet subspace t/ 3 + we present the ground state and the first three excitations, i.e. 1 < v < 4. 
The numerical results together with the available data from the literature are given in Table 4. Much more 
data are available in the literature than for the singlet states. These numbers convincingly demonstrate that 
our approach to the solution of the two-electron problem in a strong magnetic field is superior to any other 
method existing in the literature. 

Among the three related triplet states with S z = 0, ±1 the one with S z = — 1 possesses the lowest energy 
due to the spin shift BS Z . Therefore, we have given in Table 4 those lowest energies. Although for B = 
the singlet state 1 1 + is the global ground state and the triplet state l[J0 + (i.e. S z = 0) stays above it 
for any field strength, the spin shift (— B) causes the related triplet state 1 3 _ 1 + (S Z = —1) to become the 
ground state within the M = subspace above B ~ 1.112a. -a.. For the triplet ground state it is within our 
approach much easier to obtain accurate results than for the singlet ground state. This can be understood in 
the following picture: in the triplet ground state only one of the electrons occupies the one-particle ground 
state with the cusp whereas the other one occupies already predominantly an excited one-particle state and 
thus gives rise to a lower correlation contribution than in the singlet ground state. 

In analogy to the case of the singlet states we provide in Fig.l the dependence of the one-particle 
ionization energies on the field strength for the triplet states. The reference threshold T(B) has the same 
value as for the singlet case because, following our general definition of the threshold, we keep the spins 
fixed. 

We observe that the ionization energy of any state t^ 3 + in Fig.l behaves similarly to the ionization 
energy of its singlet counterpart (u + 1) 1 + . This relationship has its origin in the fact that two states 
differing only with respect to their spin symmetry possess according to the field free notation in Table 1 



similar contributions to their energies. The only difference with respect to the matrix elements of eqs.(20) 



(22) for two such states is the sign of the exchange terms. Those might be small if the particles are 
sufficiently spatially separated as it is the case for all states apart from the 1 1 + state. This statement is 
confirmed by the fact that the singlet-triplet splitting is particularly small between the 4 1 + state and the 
3 3 + state. These states belong to the field free states 3 Do and 3 3 Do (see Table 1), respectively, and the 
one-particle orbitals with d-symmetry involved here are even more spatially separated from the participating 
s-orbitals than the p-orbitals playing a role in other states with larger singlet-triplet splitting. Of course 
this argumentation breaks down when the magnetic field destroys the spherical symmetry sufficiently and 
causes a mixture of total angular momenta. Consequently, the splitting between the 4 X + state and the 
3 3 + state grows considerably at B ~ 0.02. 

4.3.2 Results for M = and odd z-parity 
a) Results for the singlet states z/ 1 0~ 
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For the singlet subspace z^O - we present the ground state and the first four excitations, i.e. 1 < v < 5. 
For B < 0.08, even the energies for the state 6 1 0~ show an excellent convergence. The corresponding data are 
presented in Table 5. For a graphical representation we choose the ionization energies which are calculated 
from the total energies with respect to the same threshold T(B) as in the case of the + subspace. The 
reason is that in a one-electron ionization process of helium constrained to keep 0~ symmetry the ionized 
electron can adopt the negative z-parity, allowing the remaining electron to occupy the one-particle ground 
state which possesses + symmetry. The ionized electron, in turn, has the same Landau energy -^ as in the 
+ case because the ^-parity does not have any influence on the Landau energy which is only assigned to 
the transversal degrees of freedom. 

The curves for the ionization energies are shown in Fig. 2. We observe that for very low fields the 6 1 _ 
state and the 5 1 0~ state are approximately degenerate which is in agreement with the fact that for two values 
off,// within the same bracket of the sequence (1), (2), (3,4), (5, 6), (7,8,9), (10, 11, 12)... the two states v 1 ^ 
and /x 1 _ correspond to the same quantum number L according to Table 1. As in the + subspace, this 
approximate degeneracy is removed for fields above B ~ 0.02. The region between B ~ 0.02 and B ~ 0.1 
exhibits some avoided crossings. 

b) Results for the triplet states u 3 0~ 

In contrast to the singlet case there exist much more data in the literature for the triplet state. In Table 
6 we have listed them together with our results for the states ^ 3 _ 1 _ (i.e. S z = —1), 1 < v < 5. Again our 
results are better than any references for finite field strengths. 

The ionization energies for the i^O" states are also shown in Fig. 2. As in the singlet case, we achieved a 
good convergence even for the fifth excited state 6 3 0~ below B = 0.08. The singlet-triplet splitting between 
the states v 1 0~ and ^ 3 0~ behaves similar as a function of the field like the corresponding splitting in the 
1 ' 3 + subspaces. 

4.4 Transitions 

In order to interpret experimental spectra from magnetic white dwarfs like GD229, it is necessary to de- 
termine transition energies from our total energies. The selection rules for electric dipole transitions are 
AS = 0, AS Z = for the spin degrees of freedom and AM = 0, AIT = ±1 or AM = ±1, AIL. = for 
the spatial degrees of freedom. With the data presented we are able to treat the AM = transitions. We 
obtain a spectrum of 30 transition energy curves for singlet transitions and 24 ones for triplet transitions 
which we show in Fig. 3a and 3b, respectively. In both logarithmic representations, we observe singularities 
belonging to zeros in the transition energies which arise due to level crossings of the initial state and the 
final state. 

Since the magnetic field of a white dwarf is not constant but varies by a factor of two for a dipole geometry, 
the spectrum of wavelengths is in general smeared out. Transitions whose wavelengths are stationary with 
respect to the magnetic field, however, reflect themselves by characteristic absorption edges in the observable 
spectrum. 

Due to this important role of the stationary lines we have summarized all transitions showing stationary 
points in Tables 7 and 8. The position and the wavelength of each stationary point have been determined 
by interpolation using the relatively crude grid of our calculations. This allows us to perform a comparison 
of the stationary components with the observed spectrum of the magnetic white dwarf GD229 |1(|. Indeed, 
the accurate data on many excited states presented here are part of an analysis accomplished very recently 
[pTCf l which clearly shows that there is strong evidence for the existence of He in the atmosphere of GD229. 
This might also be the case for other magnetic white dwarfs and therefore the present data will serve 
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astrophysicists for further comparison with observational data. 

5 Conclusions and Outlook 

We have investigated the electronic structure of the helium atom in a magnetic field by a fully correlated 
approach. We assumed the nucleus to possess infinite mass but provided scaling laws how to connect our 
fixed-nucleus results to the data which should be expected for the true finite nuclear mass. One of the goals 
of our work, the identification of the features in the spectrum of the white dwarf GD229 with electronic 
transitions in atomic helium, has already successfully been demonstrated [ 10 1. This has been possible due 



to the high accuracy of our calculations on many excited states considering the electronic structure for 
magnetic fields for < B < lOOa.u., i.e. < B < 2.3505 • 10 7 Tesla. 

The starting point of our ah initio treatment of the helium atom in a magnetic field was its full fixed- 
nucleus Hamiltonian which possesses four conserved quantities: the total spin S 2 and its z-component 
S z , the z-component L z of the electronic angular momentum and the electronic z-parity H z . Due to the 
magnetic field the spherical rotational invariance is broken resulting in the fact that the total electronic 
angular momentum L 2 fails to be a conserved quantity for nonvanishing field B ^ 0. Our full configuration 
interaction approach has been able to overcome the difficulty that the symmetry of the system changes from 
purely spherical for B = to mainly cylindrical for high fields. To this purpose we used a one-particle basis 
set of anisotropic Gaussians furnished with monomials in the coordinates p or z transversal or longitudinal 
to the field axis, respectively. The degree of anisotropy has been obtained as a result from the direct 
optimization of the nonlinear parameters for the transversal or longitudinal Gaussian by the requirement to 
solve optimally the one-particle problem of the H atom or the He + ion in a magnetic field of given strength. 
The one-particle basis functions were composed to two-particle states of pure total symmetry in the sense 
of the four conserved quantities mentioned above. In the present paper, we have provided results for the 
singlet and triplet states in the subspaces with vanishing spatial magnetic quantum number M = and 
positive or negative z-parity. 

In each of these subspaces we have built up a matrix for the Hamiltonian and for the overlap between the 
nonorthogonal two-particle states which provided a variational estimation for the energy eigenvalues after 
the diagonalization of a generalized eigenvalue problem. By very elaborate techniques for the evaluation 
of the matrix elements we were able to treat matrix dimensions of about 4000 within less than one day 
CPU time on a moderate Silicon Graphics workstation. The relative accuracy of the results for the energies 
ranged between 10~ 4 for the singlet ground states and 10~ 5 for the triplet states or general excited states 
for 2 < v < 5. To achieve this, a careful selection of the basis functions was necessary, i.e. combinations 
of one-particle magnetic quantum numbers m = ±1, ±2, ... were important to describe correlation whereas 
one-particle functions describing high excitations of the one-particle systems H or He + provided the major 
contribution also for high excitations of the two-particle system He. 

To enable a comparison of our calculated transitions with the observed spectra of magnetic white dwarfs, 
we have determined all the stationary points associated to the transitions between states with M = and 
different ^-parities. It turned out that the stationary points in the regime B ~ 0.1 and B ~ 0.3 can be 
identified with lines in the spectrum of the white dwarf GD229 which is even further confirmed by transitions 
involving M = —1 and positive z-parity[[H],|4l]]. This provides strong evidence for the existence of helium 
in the atmosphere of this white dwarf. The present data will additionally be helpful in order to see whether 
there is evidence for helium also in other magnetic cosmic objects. 

In order to complete the treatment of the electronic structure of helium in a magnetic field, we will in 
the near future calculate data for subspaces with higher magnetic quantum numbers for both odd and even 
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z-parity. Furthermore, a detailed investigation of the oscillator strengths as a function of the magnetic field 
is necessary in order to estimate properly the intensities of the transitions. 
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A The overlap matrix elements 

In the following we only give the results for the analytic representation of the matrix elements. The evaluation 
of the overlap integral is very simple because it factors into an integral over the transversal coordinate p 
and an integral over the longitudinal coordinate z. Both integrals can be reduced to a Gaussian type 

00 2 

J £™e _7i , n > 0, yielding for the dependence of the one-particle overlap (i\k) on the parameters of the 



states |i) and \k): 
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where we define (—1)!! := 1 and 7^ := 7, + j/~- The Kronecker Delta symbols reflect the orthogonality of 
two one-particle states with different magnetic quantum numbers or different z-parities. We observe that 



the overlap is real and all the parameters of the wave functions \i) and \k) enter symmetrically in eq.(24) 



B The one-particle matrix elements 

B.l The matrix elements of the kinetic energy 

The matrix elements of the operator (p + ^B x r) 2 are evaluated by using cylindric coordinates in which 
we consider 

i„ , 2 1 ( ( o 2 1 a\ id 2 d 2 ._ a 1 , ., 




(25) 

where T p ,T v ,T z represent the Laplacian, Tzeeman is the Zeeman term and T^ia the diamagnetic term. The 
evaluation of any of the matrix elements {i\Ti\k) is relatively simple since the derivatives generate prefactors 
but leave the types of the integrals unchanged in comparison to the overlap integral. A consequence is that 
the results are found to be proportional to the overlap (i\k) between the same two states. In detail we have 
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(i\T z \k) = (i\k) • { "^( 1 -"^)^ 2 +( 2 "^^+^- 1 )ftfe+"^( 1 -"^)^ . n zik ^l} (28) 

(i\T Z eeman\k) = (l\k) • -m fej B (29) 

(i\T dia \k) = (,|fc).^±l.Ii? 2 (30) 

Again we observe that each matrix element is real and symmetric with respect to an interchange of the 
states \i) and \k). Due to the selection rules in the prefactor (i\k) the occurence of the single quantity m& 
does not destroy this symmetry. Physically, the proportionality to (i\k) means that the operator of the 
kinetic energy does not couple one-particle states involving different magnetic quantum number or z-parity. 

B.2 The matrix elements of the electronic Coulomb interaction with the nucleus 

The evaluation of the matrix elements of the one-particle operator Vj = ^ is much more complicated than 
the other one-particle matrix elements discussed above. The reason is that the Coulomb potential possesses 
a spherical symmetry rather than the cylindrical symmetry of the basis functions. We have overcome 
this difficulty by applying a Singer transformation^], leaving for the spacial integrations the convenient 
Gaussian types. The remaining integration due to the Singer transformation, however, is not as simple but 
it can be solved by involving the Gaussian hyper geometric function 2-P1, yielding 

m\k) = -m- Pl ^^ rj) 2 M^ +1 > *■ ^4^+1; i-fe) (3D 

It is important to remark that very elaborate techniques are necessary to evaluate the function 2-F1 for the 
various occuring arguments with a high accuracy and with an acceptable efficiency. The standard expansion 
of 2-F\(a> b, c; z) in a power series of z is by no means sufficient because a singularity of 2F\{a,b,c;z) for 
z = 1 constrains the convergence domain of such a series to \z\ < 1. The consequence is that for ™- <C 1 

the convergence of the standard series would be arbitrarily slow. Moreover, for ™- > 2 the application of 
the standard series is completely useless. Using basis sets whose parameters ai,/3i usually cover the range 
from 10 up to 10 5 , we have thus been obliged to use various formulas for suitable analytic continuations 
of 2-F1 (O) b, c; z) to the domain \z\ > 1 ^3|. By these techniques we achieved an accuracy of 10 -13 for the 
one-particle interaction integrals with no loss of efficiency compared with the simple integrals of the kinetic 
energy. 

We remark that also the matrix elements of the electron-nucleus interaction are proportional to the 
overlap (i\k). This is in agreement with the fact that the spherically symmetric Coulomb potential Vj = - 
neither breaks the azimuthal symmetry nor disturbs the behaviour of the one-particle basis functions under 
the reflection z — > — z. 

C The two-particle matrix elements 

In contrast to the one-particle matrix elements treated above we cannot provide a single formula for the 
efficient and accurate calculation of the matrix elements {ij\Hjj\kl) of the two-particle interaction. The 
reason is the large number of different cases due to the various constellations of the individual one-particle 
functions which have to be treated differently in order to ensure a fast and highly accurate evaluation. We 
therefore provide in the following only an outline of the procedure and indicate the necessary techniques. 
Details can be obtained from the authors upon request. 
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First we apply a Singer transformation [42] in order to remove the Coulomb singularity in the integrand 
of the matrix elements, thereby introducing the new variable u according to i — - — r = —j= Jq°° du e~ u ( ri ~ r2 ) . 
Although the underlying symmetry of the basis set is cylindrical, it is advantageous to carry out the spatial 
integrations in Cartesian coordinates. The z\ ^-integral can always be factored out, and the coupling 

2 

between the two particles can be removed by the subsitution z\ — ► z\ — g u , ± . The resulting integral can be 
decomposed into a sum of [n Z ik + 1) integrals each factoring in two pure z\ and Z2 integrals of the Gaussian 
type J*^° z nz e~^ u ' z which can easily be evaluated. 

The transversal part of the electron-electron integral is much more complicated than the z-integration. 
First each term p^' |mfcl •' e -«K-mt)»'i gives rise to factors (kik + 1) ■ (\irii\ + 1) ■ (|"^fc| + 1) according 
to its decomposition to (x 2 + y 2 ) kik (x\ — i sgn(mi)yiy mi ' (x\ + i sgn(mfc)yi)l mfc l in Cartesian coordinates 
(again the same number of expressions arises for particle 2). Next, the particle decoupling transformation 
analoguous to the one mentioned above multiplies the number of integrals by an additional large factor n^. 

At this stage, all the Gaussian integrations in x\,X2,yx,y2,Zi and Z2 can be carried out, resulting in 
functions of the remaining variable u of type u 2 (l + aiu 2 ) ri (l + a2U 2 ) r2 (l + a^u 2 ) r3 {l + 04ti 2 ) r4 , where the 
fi are positive or negative integers or half-integers. The coordinates can always be chosen such that one j*j 
is a positive integer. Multiplying out the factor (1 + ajU 2 ) n enables us to reduce the remaining u-integral 
to an expression involving the Appell hypergeometric function F\(a, b, b',c; t±, £2). 



For each of the — ^ — - matrix elements (ij\Hu\kl), we thus have to evaluate n^ ■ {(fejfe + 1) • (\mi\ + 



N(N+1) 
2 

1) • (\mk\ + l)} 2 • {i^zik + 1) ' ( r i + 1) times the function F\ with in general various different arguments 
a, b, b\c. Only the two arguments t±,t2 are universal for a fixed matrix element: t\ = g +r. 1 an d £2 = 

1 — \ Q lk , J l l lk jl . For I £i| ^C 1 and |i2| <^. 1 we implemented the usual standard series expression of 



(ftfc+/3ji) Q ifc Q ji ' 



oc 



Fi(a,6,6';c;ii,i2) = E™=o fefyfSy 2^1(0 + m, 6'; c + m;£ 2 ) if, where 2 F 1 (a, 6; c; z) = E fe#fg z n is 

n=0 

the Gaussian hypergeometric function. However, the fact that for many basis functions either the argument 
|ii| or [£2] or even both happen to lie close to 1 or above makes this standard representation useless since 
its convergence domains are the unit circles with respect to £1 and £2) 44]. Therefore, it was necessary to 



use a large number of analytic continuation formulas for F\ each of which is valid for one specified class of 
arguments a,b,b',c. Even more, we have systematically compared the CPU times for altogether more than 
50 ways how to evaluate F\ and selected the fastest one for each class of arguments, involving e.g. all of the 
continuation formulas for the Gaussian hypergeometric function 2-^1 ( a ; b, c; z) given in eqs. (15.3.3-15.3. 12) 



in ref.[43]. We have additionally derived new formulas for 2-Fi(a, b, c; z) adapted for parameter constellations 
not covered in ref.|43|], and we investigated the different possibilities how to reduce F\ to 2F1 or one of its 
derivatives. 

We point out that without such a systematic analysis of the convergence properties of series represen- 
tations for F\ the present work would not have been possible: The reduction of CPU time for the same 
accuracy (10 -12 ) of our representation of the series compared to the most primitive selection of continuation 
formulas is about 10 3 . 
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Fig. 2. Ionization energies of the singlet and triplet states fc 1 and n 3 , k = 1, ..., 6, n = 1, ..., 6. The solid 
lines correspond to the singlet states, the dashed ones show the triplet states. 



J>v Fig. 3 (a) Transition energies for the singlet transitions jj, + — > rO , /i = 1, ...,5, n = 1, ...,6 as a function 
i-£h \ of the field strength. The curves are interpolated on a grid of field strengths given by our calculations. 

(b) Transition energies for the triplet transitions /U 3 + — » i^ 3 _ , \x = 1, ...,4, n = 1, ...,6 as a function of the 
field strength. 
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Table 1: Correspondence table between the field notation v +l M^ 1 ' z and the familiar field free 
notation n 2S+1 LM for the lowest 55 singlet states and the lowest 54 triplet states of He. We have 
also provided the precise total zero- field energy of each state given in the literature. 



field free 


field 


Energy 


n 2S + l L M 


^+i M (-i) u » 


-2.903724377" 


1% 


l x 0+ 


-2.145974046" 


2 1 S 


2 X 0+ 


-2.1238430864" 


2 l P Q 

2 1 P±i 


i^o- 

^(±1)+ 


-2.061271988" 


3 So 


3*0+ 


-2.0556207328" 


&D 
&D ±1 

3 X D ±2 


4*0+ 

li(±l)" 
1 1 (±2)+ 


-2.0551463620" 


3 X P 

3 x P±i 


2 X 0- 

2 1 (±1)+ 


-2.033586717" 


4% 


5*0+ 


-2.0312798461" 


A 1 D 

A l D ±l 
4 1 £>±2 


6*0+ 

2 1 (±1)- 
2 1 (±2)+ 


-2.0312551443" 


4^0 

4^+1 

4 1 F ±2 

4 1 F ±3 


3 X 0- 

3 1 (±1)+ 

iH±2)- 

1 1 (±3) + 


-2.0310696504" 


4^0 
4 1 P ±1 


4 X 0- 
4 1 (±1)+ 


-2.0211137 c 


5% 


7 X 0+ 


-2.0200158361" 


&D 

5 l D ±1 
5 1 D± 2 


8*0+ 

3 1 (±2)+ 


-2.0200029371" 


5^0 

5^+1 

5 1 i J ± 2 

5^+3 


5 X 0- 

5 1 (±1)+ 
2 1 (±2)" 
2 1 (±3)+ 


-2.0200007108" 


5^0 

&G& 
5 1 G± 4 


9 X 0+ 

4 X (±1)- 
4 X (±2)+ 

^(±3)- 

^(±4)+ 


-2.0199059899" 


5^0 
&P±i 


6 X 0- 

6 1 (±1)+ 



field free 


field 


Energy 


n 2S+1 L M 


u 2S+i M (-iF* 


— 


— 




-2.175229378" 


2 3 S 


l a 0+ 


-2.133164191" 


2 3 P 

2 3 P±i 


1 3 0~ 

1 3 (±1)+ 


-2.068689067" 


3 3 S 


2 3 + 


-2.058081084" 


3 3 P 
3 3 P±i 


2 3 0~ 
2 3 (±1)+ 


-2.055636309" 


3 3 A) 
3 3 Ati 

3 3 D± 2 


3 3 0+ 
1 3 (±1)- 
1 3 (±2)+ 


-2.036512083" 


4 a 5 


4 3 0+ 


-2.032324354" 


4 3 P 
4 3 P±i 


3 3 0~ 
3 3 (±1)+ 


-2.031288847" 


4 3 A) 

4 3 Ati 

4 3 L> ±2 


5 3 0+ 

2 3 (±1)- 
2 3 (±2)+ 


-2.031255168" 


4 3 F 
4 3 F ±1 

4 3 F ±2 
4 3 F±3 


4 3 0~ 

4 3 (±1)+ 
1 3 (±2)~ 
1 3 (±3)+ 


-2.02261624 c 


5 3 S 


6 3 0+ 


-2.020551187" 


5 3 P 
5 3 P±i 


5 3 0~ 
5 3 (±1)+ 


-2.020021027" 


5 3 A) 

5 3 Afci 

5 3 D± 2 


7 3 0+ 
3 3 (±1)- 

3 3 (±2)+ 


-2.020002957" 


5 3 F 
5 3 F ±1 
5 3 F ±2 
5 3 F±3 


6 a 0~ 

6 3 (±1)+ 

2 3 (±2)~ 
2 3 (±3)+ 


-2.020000711" 


5 3 G 

5 3 G ±1 

5 3 G ±2 

5 3 G ±3 

5 3 G ±4 


8 3 0+ 

4 3 (±1)- 
4 3 (±2)+ 
1 3 (±3)- 
1 3 (±4)+ 



" Baker et al. (1990) || 
6 Drake et al. (1992) BJ 
c Accad et al. (1971) M 



Table 2: Total energies E of the singlet ground state 
1 1 + , one-electron ionization threshold T and, if available, 
best energy values given in the literature, as a function of 
the magnetic field strength B 



B 


E(l l 0+) 


literature 


T 


0.0000 


-2.903351 


-2.903724377 a 


-2.000 


0.0008 


-2.903346 




-1.999599960 


0.004 


-2.903342 




-1.997999000 


0.008 


-2.903340 




-1.995995995 


0.020 


-2.903270 




-1.989975001 


0.040 


-2.903036 




-1.979900008 


0.080 


-2.902083 


-2.90195 d 


-1.959600176 


0.160 


-2.898290 




-1.918402804 


0.240 


-2.892035 




-1.876414090 


0.400 


-2.872501 


-2.8714 d 


-1.790105922 


0.500 


-2.855859 


-2.85385 e 


-1.734628064 


0.800 


-2.787556 


-2.77585 d 


-1.561526260 


1.000 


-2.729508 


-2.7272 e 


-1.440989741 


1.600 


-2.507952 




-1.058421519 


2.000 


-2.329780 


-2.3266 e 


-0.788842154 


5.000 


-0.574877 




1.456132354 


10.000 


3.064582 




5.609851957 


20.000 


11.267051 




14.47840453 


50.000 


38.076320 




42.45369755 


100.000 


84.918313 




90.43945348 



a Baker et al. (1990) 
d Thurner et al. (1982) 
e Larsen (1979) 
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Table 3: Total energies E of the excited singlet states v 1 + for 2 < v < 5 as a function of the 
magnetic field strength B as well as field free reference values 



B 


E(2 1 0+) 


E(3 1 0+) 


J B(4 1 0+) 


E(5 l 0+) 


0.0000 


-2.145912 


-2.061255 


-2.055613 


-2.033579 


0.0000 (lit.) 


-2.145974046 a 


-2.061271988 a 


-2.055620732852245 6 


-2.033586717026° 


0.0008 


-2.145911 


-2.061246 


-2.055607 


-2.033531 


0.004 


-2.145869 


-2.061028 


-2.055489 


-2.032867 


0.008 


-2.145743 


-2.060374 


-2.055122 


-2.031451 


0.020 


-2.144852 


-2.056838 


-2.052030 


-2.026808 


0.040 


-2.141787 


-2.050850 


-2.038691 


-2.018836 


0.080 


-2.130812 


-2.036626 


-2.004123 


-1.995610 


0.160 


-2.096743 


-1.999728 


-1.962448 


-1.945128 


0.240 


-2.054961 


-1.956971 


-1.919910 


-1.902875 


0.400 


-1.965307 


-1.865411 


-1.830581 


-1.814978 


0.500 


-1.908671 


-1.807640 


-1.773700 


-1.758191 


0.800 


-1.736422 


-1.632339 


-1.599198 


-1.584566 


1.000 


-1.617870 


-1.511771 


-1.478551 


-1.464143 


1.600 


-1.241729 


-1.130213 


-1.096276 


-1.081100 


2.000 


-0.975861 


-0.861382 


-0.826959 


-0.811912 


5.000 


1.252363 


1.379993 


1.416706 


1.432270 


10.000 


5.393200 


5.530945 


5.569418 


5.585430 


20.000 


14.248991 


14.396817 


14.437004 


14.453455 


50.000 


42.207510 


42.368689 


42.411075 


42.428159 


100.000 


90.180690 


90.351969 


90.395965 


90.414700 



a Baker et al. 
b Drake et al. 



(1989) |0 1 
(1992) |3l 



Table 4: Total energies E of the triplet states v 3 _ 1 + (S z = — 1) for 1 < v < 4 as a function of the magnetic 
field strength B. We have also provided the results given in the literature so far. 





E0 


3 0+) 


E(2 3 0+) 


£(3 3 0+) 


£(4 3 0+) 


B 


this work 


literature 


this work 


literature 


this work 


lit. 


this work 


literature 


0.0000 


-2.175220 


-2.175229 


-2.068687 


-2.068689 


-2.055630 




-2.036511 


-2.036512 






378 a 




067472° 


-2.055636309453 b 




083098° 


0.0008 


-2.176019 


-2.175050/ 


-2.069480 


-2.069278 4 


-2.056424 




-2.037278 




0.004 


-2.179190 


-2.178220/ 


-2.072505 


-2.072302 4 


-2.059507 




-2.039893 




0.008 


-2.183098 


-2.1830(5) 9 


-2.075969 


-2.0760(1)9 


-2.063153 




-2.042333 


-2.0413(4)9 


0.020 


-2.194461 




-2.084626 




-2.072654 




-2.048567 




0.040 


-2.212236 


-2.211222/ 


-2.096580 


-2.092686 4 


-2.083143 




-2.060290 




0.080 


-2.243958 


-2.2438(3) 9 


-2.121107 


-2.1209(1)9 


-2.087409 




-2.082249 


-2.0687(9)9 


0.160 


-2.296318 




-2.166290 




-2.123762 




-2.105446 




0.240 


-2.339560 




-2.206689 


-2.1927/ 


-2.162609 




-2.143909 




0.400 


-2.412723 


-2.4112/ 


-2.279523 


-2.2615/ 


-2.235815 




-2.217383 




0.500 


-2.454347 


-2.4528' 1 


-2.321985 




-2.279382 




-2.261303 




0.800 


-2.573615 


-2.5737(3)9 


-2.443352 


-2.4395(9)9 


-2.403631 




-2.386891 


-2.3497(21)9 


1.000 


-2.650655 


-2.6492^ 


-2.520902 




-2.482164 




-2.465921 




1.600 


-2.867620 


-2.8669(5)9 


-2.736687 


-2.7330(8)9 


-2.698755 




-2.682563 


-2.6665(26)9 


2.000 


-2.999708 


-2.9982^ 


-2.867135 




-2.829137 




-2.813101 




5.000 


-3.768199 


-3.7667^ 


-3.624558 




-3.584967 




-3.568584 




10.000 


-4.627450 


-4.626' 1 


-4.473459 




-4.432177 




-4.415350 




20.000 


-5.772448 


-5.771' 1 


-5.607619 




-5.564582 




-5.547290 




50.000 


-7.815256 


-7.814 h 


-7.635847 




-7.590520 




-7.563695 




100.00 


-9.843074 


-9.842' 1 


-9.652632 




-9.605634 




-9.586289 





a Baker et al. (1989) [40] 
b Drake et al. (1992) [31] 
/ Jones et al. (1996) [17] 
9 Jones et al. (1997) [22] 
h Ivanov (1994) [16] 
1 Thinner et al. (1993) [13] 



Table 5: Total energies E of the states ^ 1 0~, 1 < v < 6 as a function of the magnetic field strength B as well 
as field free reference values and the values for finite field, if available. For B > O.la.u., the energy values for 
the 6 1 0~ state are not optimally converged 



B 


E(i l o~) 


E{2 l Q-) 


E(3 l 0-) 


E(4 L Q-) 


E(5 l 0-) 


E(6 l Q-) 


0.0000 


-2.123412 


-2.055015 


-2.031252 


-2.031014 


-2.019999 


-2.019877 




-2.123843430864 6 


-2.0551463620 6 


-2.0312551443 6 


-2.03106965 6 


-2.020002937 6 


-2.019905990 6 


0.0008 


-2.123413 


-2.055009 


-2.031238 


-2.030994 


-2.019958 


-2.019824 


0.004 


-2.123388 


-2.054868 


-2.030983 


-2.030463 


-2.019427 


-2.018337 


0.008 


-2.123314 


-2.054438 


-2.030349 


-2.028842 


-2.018198 


-2.014716 


0.020 


-2.122782 


-2.051758 


-2.026727 


-2.020802 


-2.012781 


-2.005504 


0.040 


-2.120940 


-2.044641 


-2.017663 


-2.006440 


-2.000816 


-1.994681 


0.080 


-2.114310 


-2.027515 


-1.997336 


-1.983478 


-1.976059 


-1.970825 


0.160 


-2.093203 


-1.990787 


-1.957118 


-1.942289 


-1.934541 




0.240 


-2.065812 


-1.952090 


-1.916192 


-1.900745 


-1.892719 




0.400 


-2.000350 


-1.870342 


-1.831436 


-1.815130 


-1.806747 




0.500 


-1.954646 
-1.95455 6 


-1.816921 


-1.776662 


-1.759971 


-1.751519 




0.800 


-1.802966 


-1.648143 


-1.605029 


-1.587531 


-1.578773 




1.000 


-1.692794 
-1.69185 6 


-1.529609 


-1.485166 


-1.467297 


-1.458395 




1.600 


-1.331851 


-1.151043 


-1.103929 


-1.085324 


-1.076152 




2.000 


-1.072194 
-1.06905 6 


-0.883229 


-0.834929 


-0.816003 


-0.806576 




5.000 


1.137045 


1.355762 


1.408108 


1.428117 


1.437859 




10.000 


5.270657 


5.506354 


5.560836 


5.581403 


5.591347 




20.000 


14.124717 


14.372781 


14.428726 


14.449668 


14.459766 




50.000 


42.086569 


42.346227 


42.403449 


42.424712 


42.434863 




100.000 


90.064304 


90.330942 


90.388895 


90.410339 


90.420669 





6 Drake et al. (1992) [31] 
e Larsen (1979) [20] 



Table 6: Total energies E of the triplet states ^£ 1 (S z 
provided the results given in the literature so far. 



-1) for 1 < v < 6 as a function of the magnetic field strength B. We have also 





E{1 


3 o-) 


E(: 


! 3 0-) 


E(3 3 Q-) 


E(4 3 Q-) 


E(5 3 0-) 


E(6 3 0~) 


B 


this work 


literature 


this work 


literature 


this work 


literature 


this work 


It. 


this work 


It. 


this work 


It. 


0.0000 


-2.132910 


-2.133164 


-2.058016 


-2.058081 


-2.032298 


-2.032324 


-2.031252 




-2.020538 












191 6 




084 h 




354 6 


-2.031255168 b 


-2.020551187 6 


-2.020002957'' 


0.0008 


-2.133710 


-2.1323/ 


-2.058810 


-2.0584/ 


-2.033081 




-2.032038 




-2.021293 




-2.020753 




0.004 


-2.136889 


-2.1344/ 


-2.061884 


-2.0615/ 


-2.035896 




-2.034870 




-2.023839 




-2.022647 




0.008 


-2.140826 


-2.1408(8) ff 


-2.065499 


-2.0655(5)9 


-2.039077 


-2.0350(10)9 


-2.037514 




-2.026642 




-2.023039 




0.020 


-2.152378 




-2.075054 




-2.047603 




-2.041411 




-2.033564 




-2.025840 




0.040 


-2.170822 


-2.169276 4 


-2.088326 


-2.086102 4 


-2.059120 




-2.046510 




-2.041394 




-2.035472 




0.080 


-2.205130 


-2.2050(6) 9 


-2.111478 


-2.1100(3)9 


-2.079242 


-2.0578(14)9 


-2.064505 




-2.056648 




-2.051347 




0.160 


-2.266575 


-2.262806 4 


-2.154912 


-2.112992 4 


-2.118891 




-2.103193 




-2.095063 








0.240 


-2.322032 


-2.3198/ 


-2.196517 


-2.1868/ 


-2.157959 




-2.141609 




-2.133225 








0.400 


-2.422361 


-2.4196/ 


-2.275521 


-2.2463/ 


-2.233362 




-2.216039 




-2.207275 








0.500 


-2.480172 




-2.322570 




-2.278711 




-2.260925 




-2.252044 








0.800 


-2.638222 


-2.6397(9) 9 


-2.455054 


-2.4545(6)9 


-2.407425 


-2.3925(7)9 


-2.388624 




-2.379366 








1.000 


-2.733813 




-2.537231 




-2.487761 




-2.468471 




-2.459030 








1.600 


-2.987185 


-2.9873(8)9 


-2.760330 


-2.7596(8)9 


-2.706990 


-2.6749(6)9 


-2.686687 




-2.676877 








2.000 


-3.135142 




-2.893357 




-2.838224 




-2.817461 




-2.807405 








5.000 


-3.959235 




-3.657879 




-3.596162 




-3.573739 




-3.563137 








10.000 


-4.848590 




-4.509715 




-4.444111 




-4.420729 




-4.409789 








20.000 


-6.011488 




-5.645191 




-5.576755 




-5.552682 




-5.541487 








50.000 


-8.059466 




-7.672999 




-7.602384 




-7.577781 




-7.566423 








100.00 


-10.079973 




-9.688173 




-9.616899 




-9.592135 




-9.580675 









6 Drake et al. (1992) [31] 
/ Jones et al. (1996) [17] 
9 Jones et al. (1997) [22] 
1 Thurncr et al. (1993) [13] 



Table 7: Overview over all of the stationary points found in of 
the singlet transitions //0 + — » z^CT for /j, = 1,...,5, v = 1,...,5 
in the range < B < 100a. u.. The main part of the errors 
arises due to the interpolation over the relatively crude grid of 
field strengths. For high wavelengths, the finite accuracy of the 
energy values themselves contributes also to error in the energies. 



Component 


wavelength/A 


position B/a.u. 


max/min 


u 2S+1 


M (-D nz 








2 1 0+- 


-l^" 


3665 ± 13 


18.3 ±0.3 


min 


2 1 0+- 


^2^" 


4291 ± 5 


0.140 ±0.005 


min 


2 1 0+- 


^2*0" 


5169 ±11 


0.899 ± 0.003 


max 


2 1 0+- 


^3*0" 


3258 ± 5 


0.184 ±0.005 


min 


2 1 0+- 


^3*0^ 


3475 ± 6 


0.68 ± 0.02 


max 


2 1 0+- 


-^4 1 0" 


2949 ± 10 


0.17 ±0.02 


min 


2 1 0+- 


-^4 1 0~ 


3073 ± 5 


0.625 ± 0.005 


max 


2 1 0+- 


^5*0- 


2797 ± 5 


0.200 ± 0.005 


min 


2 1 0+- 


^5*0" 


2905 ± 5 


0.605 ± 0.005 


max 


3 1 0+- 


^2 l Q- 


90300 ± 3000 


0.023 ±0.001 


max 


3 1 0+- 


^2 l Q- 


43600 ± 1000 


0.11 ±0.02 


min 


3 1 0+- 


^2 l Q- 


18518 ± 100 


8.9 ±0.2 


min 


3 1 0+- 


-*&Q- 


10655 ± 40 


0.14 ±0.01 


min 


3 1 0+- 


^3*0" 


17350 ± 40 


1.45 ±0.1 


max 


3 1 0+- 


^4*0" 


7916 ± 40 


0.18 ±0.01 


min 


3 1 0+- 


^4*0" 


10250 ± 10 


1.05 ±0.10 


max 


3 1 0+- 


^h 1 ®- 


6969 ± 40 


0.18 ±0.01 


min 


3 1 0+- 


^5*0" 


8538 ± 5 


0.98 ±0.06 


max 


4 1 0+- 


^2^" 


666000 ± 190000 


0.008 ± 0.004 


min 


4 1 0+- 


^3*0" 


18000 ± 12 


0.0188 ± 0.0002 


min 


4 1 0+- 


^3*0" 


52800 ± 110 


6.6 ±0.3 


min 


4 1 0+- 


^4*0" 


13850 ± 200 


0.032 ± 0.002 


min 


4 1 0+- 


-^4 1 0" 


41630 ± 220 


1.75 ±0.01 


max 


4 1 0+- 


^5*0- 


11530 ± 100 


0.025 ± 0.002 


min 


4 1 0+- 


^5 l 0- 


22750 ± 30 


1.28 ±0.02 


max 


4 1 0+- 


-^o- 


9730 ± 12 


0.0239 ± 0.0003 


min 


5 1 0+- 


^4*0" 


25000 ± 2000 


0.059 ±0.001 


min 


5 1 0+- 


^5*0- 


34390 ± 40 


0.0085 ± 0.0015 


max 


5 1 0+- 


-t&Q- 


23000 ± 600 


0.067 ± 0.005 


min 



Table 8: Same as Table |7| but for the triplet transitions ^/ 3 H 
z^ 3 0~ for [j, = 1, 



,4, v 



1,...,6 



Component 


wavelength/A 


position B/a.u. 


max/min 


V 2S+1 


M(" 1 ) nz 








1 3 0+- 


-► i 3 o~ 


1869 ± 5 


40 ±2 


min 


1 3 0+- 


^2 3 0~ 


3184 ± 3 


0.224 ± 0.005 


min 


1 3 0+- 


^2 3 0" 


4289 ± 3.5 


2.25 ±0.02 


max 


1 3 0+- 


^3 3 0" 


2503 ± 3 


0.281 ± 0.005 


min 


1 3 0+- 


^3 3 0~ 


2837 ± 2 


1.51 ±0.01 


max 


1 3 0+- 


^4 3 0~ 


2292 ± 3 


0.293 ± 0.005 


min 


1 3 0+- 


^4 3 0" 


2522 ± 1 


1.39 ±0.01 


max 


1 3 0+- 


^5 3 0" 


2197 ±5 


0.299 ± 0.005 


min 


1 3 0+- 


^5 3 0~ 


2394 ± 1 


1.34 ± 0.02 


max 


2 3 0+ - 


^2 3 0~ 


56500 ± 2100 


0.05 ±0.01 


max 


2 3 0+ - 


^2 3 0" 


39900 ± 300 


0.17 ±0.01 


min 


2 3 0+ - 


^2 3 (T 


12090 ± 25 


26.6 ±0.3 


min 


2 3 0+ - 


^3 3 0~ 


9347 ± 3 


0.24 ±0.01 


min 


2 3 0+ - 


^3 3 0~ 


16130 ±45 


3.65 ±0.15 


max 


2 3 0+ - 


^4 3 0~ 


6990 ± 10 


0.26 ±0.01 


min 


2 3 + - 


^4 3 0" 


9180 ± 5 


2.23 ± 0.06 


max 


2 3 0+ - 


^5 3 0~ 


6186 ± 10 


0.27 ±0.01 


min 


2 3 0+ - 


^5 3 0~ 


7630 ± 3.5 


1.87 ±0.05 


max 


3 3 0+- 


^3 3 0~ 


18000 ± 2000 


0.026 ± 0.003 


min 


3 3 0+- 


^3 3 0~ 


37428 ± 80 


20.5 ±0.5 


min 


3 3 0+- 


^4 3 0" 


12400 ± 1000 


0.039 ± 0.003 


min 


3 3 0+- 


^4 3 0~ 


22200 ± 100 


0.14 ±0.02 


max 


3 3 0+- 


^4 3 0~ 


21650 ± 50 


0.25 ±0.02 


min 


3 3 0+- 


^4 3 0" 


40600 ± 40 


4.5 ±0.3 


max 


3 3 0+ - 


^5 3 0~ 


10400 ± 300 


0.030 ± 0.002 


min 


3 3 0+- 


^5 3 0" 


15950 ± 100 


0.135 ±0.015 


max 


3 3 0+- 


^5 3 0~ 


15450 ± 50 


0.275 ±0.015 


min 


3 3 0+- 


^5 3 0~ 


21000 ± 55 


2.6 ±0.3 


max 


3 3 0+ - 


^6 3 0" 


9400 ± 100 


0.032 ± 0.002 


min 


4 3 0+ - 


^3 3 0~ 


477000 ± 15000 


0.0225 ±0.0015 


max 


4 3 0+ - 


^3 3 0~ 


150000 ± 90000 


0.08 ± 0.02 


min 


4 3 0+ - 


^4 3 0" 


96500 ± 1000 


0.0062 ± 0.0002 


max 


4 3 0+ - 


^4 3 0" 


20150 ± 300 


0.064 ± 0.001 


min 


4 3 0+ - 


^4 3 0" 


206000 ± 5000 


0.17 ±0.03 


max 


4 3 0+ - 


^4 3 0" 


181000 ± 12000 


0.26 ± 0.04 


min 


4 3 0+ - 


^5 3 0~ 


31400 ± 400 


0.014 ±0.001 


max 


4 3 0+ - 


^5 3 0~ 


16800 ± 500 


0.068 ± 0.005 


min 


4 3 0+ - 


^5 3 0" 


44700 ± 3000 


0.18 ±0.02 


max 


4 3 0+ - 


-*5 3 0~ 


38000 ± 6000 


0.31 ±0.03 


min 
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